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Abstract 

We present an improved measurement of the rf meson energy spectrum in T(IS') decays, using 
1.2 fb _1 of data taken at the T(IS') center-of-mass energy with the CLEO III detector. We 
compare our results with models of rf gluonic form factor that have been suggested to explain the 
unexpectedly large B — > rfX s rate. Models based on perturbative QCD fail to fit the data for large 
rf energies, showing that Standard Model strong interaction dynamics is not likely to provide an 
explanation for the large rate of high energy rf observed in B decays. 



* Submitted to the 33 rd International Conference on High Energy Physics, July 26 - August 2, 2006, Moscow 
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I. INTRODUCTION 



A surprisingly large rate for B — > rj'X s , with high momentum rj', (p v /=2 - 2.7 GeV/c), 
was observed by CLEO 0, 0] and confirmed by BaBar 3]. This result motivated intense 
theoretical activity because new physics could account for such an enhancement. However, 
Standard Model explanations have also been proposed. For example, Atwood and Soni |3] 
argued that the observed excess is due to an enhanced b — > sg penguin diagram, comple- 
mented by a strong rj'gg* coupling, induced by the gluonic content of the rf wave function. 
Fig. ^ (left) shows the corresponding Feynman diagram. The high q 2 region of the q*gv' 
vertex function involved in this process also affects fast rj' production in T(IS') decay p, |6(, 
whose relevant diagram is shown in Figure ^ (right). Thus a precise measurement of the 
rj' inclusive spectra from the process T(IS') — > ggg* — > rj'X provides valuable information 
towards our understanding of important B meson decays. 
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FIG. 1: Feynman Diagram for b — » sg (left) and T(1S) — > ggg* — » rj'X (right). 

The effective vertex function r/'gg* can be written as 

H(q 2 )6 ab e a ^q a k e1e^ 



(1) 



where q is the (g*) virtual gluon's four-momentum, k is the (g) "on-shell" gluon's momentum 
(k 2 = 0), a, b represent color indices, £i, e£ are ^ ne polarization vectors of the two gluons, 
and H(q 2 ) is the g*gr/ transition form factor. Different assumptions on its q 2 dependence 
have been proposed 0, 0, H i, E3 • 

While ARGUS was the first experiment to study the inclusive rj' production at the T(IS') 
HH, cleo ii H3 was the first experiment to have sufficient statistics to measure inclusive rj' 
production from the subprocess T(1S) — > ggo- These data ruled out a class of form factors 
characterized by a very weak q 2 dependence 0, 0| • Intense theoretical activity has followed 
0| to derive the perturbative QCD form factors from models of the rj' wave function. 
Attempts to use CLEO II data to constrain the model parameters |l3| were not conclusive, 
due to the limited statistics at the end point of the r/ spectrum. Moreover, the order of 
the perturbative expansion necessary to achieve a good representation of the data was not 
clearly defined |l3[. Thus an improved measurement, based on a higher statistics sample, 
is important to settle these issues. This work reports a new measurement of the inclusive 
T]' spectrum from the process T(IS') — > ggg* —> rj'X based on the largest T(1S) sample 
presently available, more than a factor of 11 higher than the previous study [l2j. Thus the 
measured high energy distribution function provides a much more stringent constraint. 
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II. DATA SAMPLE AND ANALYSIS METHOD 



We use 1.2 fb -1 of CLEO III data recorded at the T(IS') resonance, at 9.46 GeV center- 
of-mass energy, containing 21.2 x 10 6 events and off-resonance continuum data collected at 
center-of-mass energies of 10.54 GeV (2.3 fb -1 ). 

The CLEO III detector includes a high resolution tracking system ^ij, a state of the art 
Csl electromagnetic calorimeter |15J], and a Ring Imaging Cherenkov (RICH) hadron iden- 
tification system The Csl calorimeter measures the photon energies with a resolution 



of 2.2% at E = 1 GeV and 5% at £'=100 MeV. The tracking system provides also charged 
particle discrimination, through the measurement of the specific ionization dE/dx. 

We detect 7/ mesons through the channel 7/ — > t]tt + tt~ and 77 — > 77. The branching 
fractions for these processes are (44.5±1.4)% and (39.38±0.26)% respectively. We identify 
single photons based on their shower shape. The photon four-vectors are constrained to have 
invariant mass equal to the nominal 77 mass. Subsequently, 7/ candidates are combined with 
two oppositely charged tracks to form an 77'. Loose ir consistency criteria based on dE/dx 
measurements are applied to the charged tracks. 

The gluonic 77' production at the T(IS') is expected to be dominant only at very high q 2 , 
or, equivalently, at high 7/ scaled energy Z, where Z is defined as 

<7 _ E n > _ Erf 



Ebcam 2Mx(lS) 

where E v r is the 7/ energy and -Ebeam is the beam energy. Enhanced 7/ production at high 
Z would indicate a large r] f g*g coupling. 

For low energy t/s, photons coming from low energy 7r°s are a severe source of background. 
Thus a 7T° veto is applied for Z < 0.5, whereby photon pairs that have an invariant mass 
consistent (within 2.5 cr) of the nominal 7r° mass are not included as the candidate photons 
for 77 reconstruction. We consider only 77' with scaled energy Z between 0.2 and 1 and divide 
this range into eight equal bins. Fig. |21 demonstrates the extraction of the 7/' yields in these 
bins for the Y(IS') sample. Fig.[S]shows the corresponding plots from the continuum sample 
taken at a center-of-mass energy of 10.54 GeV. In order to derive the rj' signal yields, we 
fit the mass difference spectra AM V ' V = M(7i + tt~i]) — M(i]) in each Z bin with a Gaussian 
function representing the signal, and a polynomial background. The Gaussian is used only 
to define a ±2.5<r signal interval. Then the 7/ yield in this interval is evaluated counting 
events in the signal window, after subtracting the background estimate obtained from the 
fit function. As the signal is not described well by a single Gaussian function, this procedure 
minimizes systematic uncertainties associated with the choice of an alternative signal shape. 

Information on the gluon coupling of the 7/ can be drawn only from the decay chain 
Y(IS') — > g*gg — > rj'X. Thus we need to subtract both continuum 7/ production and 7/ from 
the process T(IS') — > 7* — ► qq. The latter component is estimated using 

B(T{IS) ^qq) = R- B{T{1S) -» = (8-83 ± 0.25)%, (3) 

where R is the ratio between the hadronic cross section 7* — ► qq and the di-muon cross 
section 7* — > at an energy close to 9.46 GeV. We use R = 3.56 ± 0.07 [ljj and 

B(T(1S) -> fi + fi-) = (2.48 ± 0.05)% |3. 

The two dominant components of the 7/ spectrum have different topologies: T(1S) — ► ggg 
produces a spherical event topology, whereas qq processes are more jet-like. This difference 
affects the corresponding reconstruction efficiencies. Fig. 0] shows the efficiencies obtained 
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FIG. 2: The difference of the rjir + TT~ and r/ invariant masses spectra in different Z ranges re- 
constructed from T(IS') data, fit with a single Gaussian function for the signal and a first order 
polynomial for the background. 

for the two event topologies with CLEO III Monte Carlo studies. We use GEANT-based 
[l9l | Monte Carlo samples, including T(1S) and continuum samples. In order to use the 
continuum sample taken at 10.54 GeV center-of-mass energy for background subtraction, 
we develop a "mapping function" , to correct for the difference in phase space and Z range 
spanned in the two samples. The procedure is described in detail in Ref . |l2| . By comparing 
the rj' energy distribution functions for the Monte Carlo continuum samples at center-of-mass 
energies equal to 9.46 and 10.54 GeV, we obtain the mapping: 

Z 9 46 = -0.215 x 10" 2 + 1.2238 Z 10 54 - 0.6879 (Z 10 54 ) 2 + 0.8277 (Z 10 54 ) 3 - 0.3606 (Z 10 54 ) 4 , 

(4) 
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FIG. 3: The difference of the r]Tr + ir~ and r/ invariant masses spectra in different Z ranges recon- 
structed from continuum data taken at a center-of-mass energy of 10.54 GeV, fit with a Single 
Gaussian function for the signal and a first order polynomial for the background. 



where Z 9AG is the Z value used to subtract the continuum contribution, as mapped from 
Z 10 .54, the measured Z in the continuum data taken at the center-of-mass energy equal to 
10.54 GeV. 

The 'jgg/ggg cross section ratio is only about 3%, thus we make no attempt to subtract 
this term. 

Upon correcting the raw yields with the corresponding efficiencies, the rf spectrum from 
the 3 gluon decay of the T(IS') is extracted using the relationship: 

N(T(1S) - ggg) = iV had - iV( 7 * - qq) - N(T(1S) - qq), (5) 

where iVhad is the number of hadronic events in our sample, and iV(7* — > qq) is the number of 
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continuum events estimated from the sample taken at 10.54 GeV, corrected for the luminosity 
difference between the resonance and continuum data, and the s dependence of the cross 
section for the process 7* — > gg. The component N(T(1S)) —> qq is estimated from 7V(7* — > 
gg) using the relationship: 

N{T(1S) - gg) = iV( 7 * -> gg) x ^isr-off/isr-on x 

= iV(7* -> gg) x .Risr-off/isr-on x , (6) 

where -Risr-off/isr-on accounts for the difference between the T(1S) — >■ gg — > rj'X and the 
7* — > gg — > spectra due to initial state radiation (ISR) effects, estimated using Monte 
Carlo continuum samples with and without ISR simulation, and a(T(lS) — > qq)/a(e + e~ — > 
qq) accounts for the relative cross section of these two processes. The correction factor 
-Risr-off/isr-on differs from 1 by a few percent at low Z and as much as 25% at the end 
; joint of the 7/ energy. The cross sections used are a(T(lS) — > H + = 0.502 ± 0.010 nb 

18^ and a(e + e~ — > ji + fi~) is 1.38 nb |2flj . Fig. EJshows the differential cross sections da/dZ 

'or the processes T(IS') — > ggg, T(IS') — > qq, and 7* — > qq. 
Theoretical predictions give the energy distribution function 

dn 1 dN(T(lS)^ggg) , 

dZ N{T(1S) -> ggg) X dZ 

we obtain the corresponding experimental quantity by dividing by the total number of 
N(T(1S) — > ggg), estimated by applying Eq. |S] without any Z restriction. Fig. EJa) shows 
the T(1S) — > ggg — > r/X energy distribution function, whereas Fig. EJb) and c) show the 
corresponding distributions for T(IS') — > qq — > rj'X, normalized with respect of the total 
number of T(IS') — > qq and Y(IS') — > t/X, normalized with respect to the total number of 
T{1S). 

The inclusive rj' production at the T(15) is expected to be dominated by the transition 
T(liS') — > ggg* — * rj'X only at high rj' energy. The energy at which this occurs cannot be 
predicted from first principle: an empirical criterion is the y 2 of the theory fit to the data. 
For example, a numerical analysis of the CLEO II data [13j obtained a x 2 of 2.4 for three 
degrees of freedom, using the 3 experimental points at Z > 0.7, and ss 24 for 4 degrees of 
freedom using the 4 points at Z > 0.6. Thus Ali and Parkhomenko concluded that the Z 
region likely to be dominated by the T(15) — > ggg* — ► rj'X starts at Z > 0.7. Thus, we 
quote separately global branching fractions for T(15) — > rj'X and the corresponding results 
for Z > 0.7. 

Table summarizes the dominant components of the systematic uncertainties. The sys- 
tematic errors on the branching fractions from rf are ±8.1% for gg, ±9.1% for the ggg — > rj'X 
for Z > 0.7 and ±7.2% for all other branching fractions. 

Thus we obtain: 

n(T(15) -+ (ggg) -+ rfX) S ^(15) = (3-2 ± 0.2 ± 0.2)%, 



n(T(15) - (gg-) - r/Jf) = = ( 3 " 8 ± °' 2 ± °- 3 ) % ' 

n(T(15) - rj'X) = = (3-0 ± 0.2 ± 0.2)%. (7) 
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Sources 


ggg Sample (Z > 0.7) 


qq Sample 


All others 


Reconstruction efficiency of tx^ 


2 


2 


2 


Reconstruction efficiency of rj 


5 


5 


5 


Number of rf from fit 


1 


1 


1 


Total number of T(1S) 


1 


1 


1 


B(n' — > ir + ir~r]) 


3.4 


3.4 


3.4 


B(T(15) - gg) 


- 


3 


- 


Ratio of integrated luminosity 


1.9 


1 




°T(l,Sr)->j4+fi- ) cr e+e--»|* + M _ 


0.7 


1.3 




7T° veto 




1.7 


0.4 


Z mapping 


6 


3 


3 


Total 


9.1 


8.1 


7.2 



TABLE I: The components of the systematic errors (%) affecting the branching fractions reported 
in this paper. 



The T(1S) — > rj'X branching fractions at high momentum (Z > 0.7) are measured to be: 

n(T(15) -> (ggg) -> t/X)) z>0 . 7 = 3.7 ± 0.5 ± 0.3) x 10" 4 , 
n(T(lS) -> (gg) ?/X)) z>0 . 7 = (22.5 ± 1.2 ± 1.8) x 10~ 4 , 

7i(T(15) -» v'X)z>a.y = (5.1 ± 0.4 ± 0.4) x 10" 4 . (8) 

III. COMPARISON WITH THEORY AND CONCLUSIONS 

A. Kagan used the ratio Rz>o.7 defined as 

(9) 

Z>0.7 

to obtain a first rough discrimination between radically different q 2 dependence between the 
form factors. At the time that Ref. [H[ was published, the comparison was based on 90% 
c.l. upper limits on the data. This test repeated with our present data give values of Rz>o.7 
equal to 74 for a representative slowly falling form factor 7], 25 for the intermediate form 
factor studied by Ref. and 2 for the perturbative QCD inspired shape. Thus the last 
shape is the closest to the data, although not providing a perfect match. 

Different approaches have been taken to implement perturbative QCD calculations. Their 
difference is translated into different assumptions for the form factor H(q 2 ), Kagan and 
Petrov |H assume H(q 2 ) fa const w 1.7 GeV -1 , Ali and Parkhomenko related H(q 2 ) to the 
expansion of the two light-cone distribution amplitudes (LCDA) describing the quark and 
gluon components of the rf wave function. Fig. [7| shows the measured dn/dZ distribution, 
compared with three representative choices for H(q 2 ). H(q 2 ) — H — 1.7GeV" _1 , H as , based 
on the asymptotic form of the rj' meson LCDA, and H(q 2 ) corresponding to the spectrum 
with the Gegenbauer coefficients from the best fit range obtained by Ali and Parkhomenko 
[l3| , using the previous CLEO II data and the constraints from the rj' — 7 transitions ■ 
Note that most of the discrepancy between theory and data occurs in the Z = 0.7 bin. In 
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Z>0.7 



n t h 
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fact, the x 2 f° r the fit of the new data with this theoretical parametrization is 27. This may 
imply that higher order terms in the QCD expansion need to be taken into account, or that 
the T(1S) — * 3g is not the dominant source of r/, at least at a scaled energy as high as 
Z = 0.7. 

In conclusion we have measured the energy spectra of the rj meson in the decay T(IS') — > 
rj'X. Our results are not very well described by existing models based on strong gluonic 
coupling of the rj'. Thus the observed B —> rj'X inclusive branching fraction is unlikely to 
be explained by an enhanced g*gr]' form factor, and an explanation outside the realm of the 
Standard Model may be needed to account for this large rate. 
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FIG. 4: The rj reconstruction efficiencies as function of Z for different MC samples with no ir c 
veto, and with 7r° veto in photon selection. The 7r° veto was applied in this analysis for Z < 0.5. 
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a) Via ggg 
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FIG. 6: The energy distribution function dn/dZ as denned in context for a) T(1S*) — > ggg — > rj'X, 
b) T(15) -» qq -> r/'X, and c)T(15) -> r/'X. 
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FIG. 7: Energy spectrum of the r/-meson in the decay T(1S) —> rj'X: measured spectra: open 
squares correspond to the data presented in this paper, filled circles are the previously reported 
CLEO II data ^3|:a) measured spectra b) Z > 0.7 region zoomed in to show the comparison 
with the theoretical predictions more clearly. The dashed curve corresponds to a constant value 
of the function H(p 2 ) = Hq ~ \.7GeV~ l , and the dash-dotted curve (H as (p 2 )) corresponds to the 
asymptotic from of the r/-meson LCDA (i-e., = and B^ = 0). The spectrum with the 
Gegenbauer coefficients in the combined best-fit range of these parameters is shown in the shaded 
region. The solid curve corresponds to the best fit values of the Gegenbauer coefficients from the 
analysis of the T(15) — > rj'X CLEO II data alone. The dotted points represent the CLEO II data 
and the boxed points represent the CLEO^! data. 



